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Abstract — The fatigue behavior of asphalt mixtures is the 
basis of the design of road structures. The heterogeneous nature 
of the material and the thermo-viscoelastic behavior of the 
bitumen could make the interpretation complex. Damage is 
defined based on the concept of effective stress associated with 
the equivalence deformation assumption. This concept translates 
both the direct and semi-direct methods. Modeling the 
mechanical behavior leads to define the materials in order to 
know its viscoelastic behavior in the frequency domain and thus, 
the complex modulus test. 

Index Terms — Fatigue, damage, thermomechanical, asphalt, 
bituminous, viscoelastic 

I. Introduction 

The bituminous coated materials are reconstituted from 
granules (0-6 mm for the studied mixture) that form the 
skeleton and whose cohesion is ensured by bitumen. The 
resulting mixture is a heterogeneous material from which the 
matrix imparts viscoelastic properties, but its characteristics 
are highly dependent on temperature. The modeling of 
characteristics’ thermal and time dependence is acquired in 
the field of small perturbations. 

A. Loading 

Although attempts have been carried out to simulate in the 
laboratory the realistic and random loading conditions [1], 
most of the time, the applied loading cycles are periodic and 
do not take into account the amplitude or frequency variations 
on real roadways. The main forms of stress signals used for 
fatigue testing are listed in Fig. 1 . 


materials, others consider it is very important to analyze the 
results of stress tests mainly through the influence of the rate 
of application of the load on the complex modulus [3]. 

B. Mechanical Properties 

The mechanical properties can be deduced from the analysis 
of results of laboratory tests. These tests, which are standard 
when used in design, are made under defined conditions of 
temperature and load. They are performed on samples made 
up in the laboratory or eventually 

The tests that characterize the mechanical properties of the 
coatings are many and varied, it is nevertheless possible to 
classify them into different categories. A first classification 
into three categories, which are more related to the 
interpretation of the test rather than its nature, is suggested by 
BONNOT (1973, 1984) [4] and [5]. A second classification 
proposed by DI BENEDETTO (1990) introduced two test 
classes: homogeneous tests and non-homogeneous tests [6]. 

C. Fatigue of bituminous mixtures 

The fatigue phenomenon of asphalt mixtures is characterized 
by its breaking after repeated application of a large number of 
loading (in practice greater than 10 4 ) whose amplitude is less 
than the resistance to the instantaneous breaking of the 
material. 

By definition, a test tube is called tired when his unit reached 
half of its initial value measured at the first loading cycle in 
the same test conditions (temperature and frequency). Life N f 
of a test tube is the number of cycles corresponding to a 
modulus of rigidity equal to half the initial modulus of the 
same test tube (Figure 1.2). 
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Figure 1: Behavior of Applied cycles used for fatigue 
testing [2] 



Figure 2: Definition of traditional criterion for determining 
the fatigue life N f 


While some authors consider that the shape of the applied 
cyclic loading has little influence on the service life of 
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The classic test is one of the most widespread fatigue criteria. 
The life Nf determined from this test is used for the design of 
pavements. 

According to the solicitation level and the number of applied 
cycles, the asphalt has three types of behavior (Fig. 3): 
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Figure 3: Tests on coated homogeneous (temperature 20 ° C) 

[7] 


- For low number of loads and deformations of a few 

percent, the test is a deformability test and a highly 
nonlinear behavior is observed. 

- For loads consisting of a few hundred cycles and low 

distortion (< 1G _ “) the behavior is considered a first 

approximation, linear viscoelastic behavior. 

- When loading several thousand cycle and under low 

deformation, damage phenomena appear, the material 
fatigue. 

Bituminous mixes are submitted, on the road, to short-term 
stresses at each passage of a heavy vehicle. Thus, the floor, 
which behaves like a rigid material resting on a deformable 
support, is subjected to bending forces. Tensile stresses then 
appear at the base of the floor and are repeated at each wheel 
arch. Since the asphalt is less resistant to traction and 
compression, in the case of passing a large number of 
vehicles, cracks are initiated on the basis of the floor due to 
the fatigue phenomen 


1 ) Timing Factors 

The reality differs from the theory and requires timing 
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Where 

• Calibration temperature , 

E 1 ^ = cst ; 

• k c correlation coefficient lab / field k c E [l; 1,3] ; 

• kp risk coefficient (reliability); 

• k s default coefficient of lift of the 

platformk E E [l; 1,2] ; 

• k^ coefficient of discontinuity and thermal gradient. 
D. Damage 

Damage as it was defined, based on the concept of effective 
stress associated with the deformation equivalence 
hypothesis, is not directly accessible to measurements, 
possible evaluation methods being [6] : 

- The semi-called direct method where we measure the 
variation of the surface or the volume of microvoids or 
density of the material; 

- Indirect methods where damage is measured through its 
effect on the coupling deformation / damage. 


1 ) Semi-direct methods 

This method requires physical observation tools 
(microscope,...) for measuring the surface or the volume of 
microvoids, and in all cases, it poses serious practical 
difficulties. Also, the theoretical difficulties appear, 
(mechanical interpretation of delicate measurements) [9] . 

2 ) Indirect Methods 

Measurements by the coupling deformation-damage resulting 
from the modification of the mechanical properties caused by 
the damage (measuring characteristics of elasticity, 
viscoelasticity, plasticity, Viscoplasticity, etc.). 

These measures are thus closely linked to the chosen model 
to describe the behavior of the material. We find in [9], [10], 
measurement techniques for elasticity, plasticity and 
Viscoplasticity. 

We intend to analyze the possibilities of measuring action 
through the viscoelastic model as it was conditioned above. 
The damage is linked to the history of the stress which makes 
complex determination in every load case. We are interested 
in sinusoidal loads which we will use later. 

- /A E(t - (t) dx (3) 

Suppose a sinusoidal deformation in stress for different stages 
of damage, it will be assumed that the stress response is 
sinusoidal (deformation fatigue test) and shifted with respect 

to deformation of an angle cp. 

So we have: 

fi(t3 =£* e CE^) = * e fc + Ct)I (4) 

with E + (t) = £. 

and a(t) = a. ^ E (e ] Llit " = J? e [a + (t)](5) 
with 

Equation (5) then gives: 

“ re ] £ C t - O rfT (6) 

This expression (7) becomes: 


We recognize there the LAPL ACE-CARS ON transform E(t) 
where the argument is ico (complex modulus of undamaged 
material). 

jft) 






Assuming that the phase difference is very little influenced by 
the phenomenon of micro cracking (assumption implicitly 
assumed in the stress tests), 1 - D(t) is considered 
non-periodic variable but slowly over time. We can write 
complex variable between: 

fF + (0 = [1 -D(t)lEXitoU + = [1 -D(t) ]-£Cw)-fl^-£ + CtX9) 


Let 

eGu,t) = [l — D(t3].EGo](i0) 

and 

E*Gto,t) = Etcui, t). e 3ip (11) 

E L represents the complex modulus of the material at time t. 
Then we have 

d w =1 -to (12) 
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The damage is to be deducted from the modification of the 
complex modulus of the material assuming the argument cp to 
be constant. This expression is only valid in isothermal 
conditions. Temperature rises influencing the complex 
modulus are to be considered in any assessment of the damage 
through this characteristic in the asphalt-coated materials. 

II. Mechanical modeling of the behavior of asphalt 

MIXTURES 

Mechanical modeling of the behavior of the material must be 
performed in view of the structure calculation. It does not 
represent the behavior of the material but its behavior in the 
structure, that’s in its context of use; loading, speed loads, 
temperatures and environment. 

A. Materials 


B. Characterization of viscoelastic behavior of asphalt 
mixtures 

The asphalt is a viscoelastic material that is to say, its 
mechanical behavior varies with the conditions of stress: 
speed, load time, temperature. To characterize the changing 
behavior of asphalt mix from the state "liquid" to that of a 
viscoelastic solid, the test measuring their complex modules 
was developed by Huet[ll] It determines the intrinsic 
characteristics of materials (modulus and phase angle) for a 
large number of temperature / frequency pairs. 

After describing this test and from standard results, 
we will see that the asphalt mixtures verify the time / 
temperature equivalence principle and it is then possible to 
know the properties of the coated at any temperature by 
building its main curve. 


1 ) Aggregates: Aspect of the grading curve 

It seems that the grading curve is not an easy parameter to 
study since it can’t be considered solely insofar as when the 
filler content and the bitumen content are a must. However, all 
the results show that continuous grading curve is better than a 
discontinuous curve. 

The European standard defines the aggregate as the "granular 
material used in construction. An aggregate can be natural, 
artificial or recycled." 

• the natural aggregate is the aggregate of mineral origin 
which has not undergone any processing other than 
mechanical. In this category we settle rock aggregates, such as 
limestone, porphyry, trap ...; 

• the artificial aggregate is the aggregate of mineral origin 
resulting from an industrial process comprising thermal or 
other changes. In this category we arrange the processed 
aggregates, such as expanded shale, expanded clay, expanded 
mica (vermiculite) ... 

• the recycled granulate is a granulate resulting from the 
processing of inorganic materials previously used in 
construction. In this category fall the aggregates such as 
crushed concrete, the milled asphalt mix... 

2) II. 1.2 Fines 

The fines, also called fines or fines addition, is a fine 
aggregate, of approximately 0-125 pm and intended to 
perform, to load various products of the construction industry 
BPW (Building and Public Works). 

It support: 

• the artificial asphalt in doses close to 10%, for its power 
stiffening; 

• bituminous coated in neighboring doses of 5 to 6% in the 
countries of northern Europe, 3-4% in the countries of 
southern Europe, for its power stiffening; 

• hydraulic cement concretes, in doses close to 50 to 80 kg per 
m 3 , as a substitution of the cement. 

3) II. 1.3 Hydrocarbon binder 

A hydrocarbon binder is generally an adhesive material (a 
binder) comprising bitumen, tar, or both. This combined 
element together with aggregates provides "embedded 
materials." 


1 ) characterization in the frequency domain 
The characterization of the behavior of bitumen in the 
frequency domain is done through the trial known as the 
complex modulus. 


That’s to apply to an asphalt specimen 0 for a given 
temperature, a sinusoidal stress of a frequency f = co / 2 tt in the 
field of small deformations. In the case of a linear viscoelastic 
behavior, Mandel [12] showed that the response to such 
request is also a sine wave but phase-shifted by an angle §. 

In practice, a deformation E"(t, to, 0„J = E c 5in([jjtJ is applied 
to a test piece of bitumen and response stress 


o m ( t, u, 0) = a 0 sin(tdt -f <f0 is recorded 
The complex modulus E"(td, 0) is defined by the ratio 
between the stress cr“(t, qj, OJand the deformation E“(t, to, 0j . 
We can write: 



(f.il. 9 ) 


(13) 


This complex module can also be written in the forms of 
equations: 

E = | 0)( g (14) 

and 

) = E L (t li.0) + EjEjCw, (15) 

Where : 

- | E * (eo) | and tfi (oo) are the norm and the phase angle of 

the complex modulus of the material under the conditions 
of the test [14], [15] ; 

- E 1 (to r 0j represents the dynamic elastic modulus which 

reflects the instantaneous response of the material. This is 
the real part of the complex modulus. It quantifies the 
stored elastic energy; Ei is reversible module associated 
with the elastic portion of the material; 

- E 2 (ifi? 0)represents the loss modulus which represents the 

viscosity (irreversible module). It quantifies the energy 
dissipated by internal friction as a result of a solicitation. 
This dissipated energy is transformed into heat which 
increase the temperature within the material under cyclic 
loading. 

2 ) trapezoidal Specimens 

As all flexural tests, it is a non-homogeneous test; the 
trapezoidal 

specimen is embedded at its large base and biased to its top. 
Solicitations may be exercised in strength or displacement. 
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The trapezoidal shape of the test piece is chosen to obtain a 
maximum deformation outside the mounting region of the 
sample loaded in cantilever beam. Rupture is generally 
performed in the vicinity of 1/5 of the total height of the 
specimen, great experimental disparity exists in the values of 
this height. 

These trapezoidal specimens have a large base of 50 mm, a 
small base of 25 mm, a height of 250 mm and a thickness of 
25 mm. [15]. This sample dimensions allow considering the 
asphalt samples as continuous and homogeneous media. 



Figure 4: (a) Trapezoidal specimen for 2PB and (b) 2 Point 

Bending Test 


The applied and the recorded waveform look symmetrical 
around the point zero. Fig. 7 shows the typical trend of the 
stiffness modulus. It is calculated considering the recorded 
stress divided by the applied strain. As it can be noticed, a 
three stage evolution process is recorded during a fatigue test. 
After a rapid evolution of stiffness (phase 1), due to the 
internal heating phenomenon, stiffness decrease seems more 
regular (phase 2). Fracture occurs in the final stage (phase 3) 
and it is characterized by an acceleration of stiffness drop 
[16]. 
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Figure 6 Stiffness evolution during a fatigue test (160pe) 



The maximum values of stress and deformation and 
E miK achieved by the specimen are calculated from geometric 
characteristics of samples using formulas from a strength of 
materials calculation. 




£|r “ r - 


,rf 


(16) 




aft 


2 aa[b-£) 




(17) 


The development of these expressions to calculate the module 
E Then allowing the calculation of E^nd E 2 . 


E i _=}'(^cos<p+v*} 1t y\8) E 1 =r(^sin4>) (19) 


Where F is the force measured, co the angular frequency, y is 
the form factor depending on the dimensions of the specimen, 
D movement at the top of the specimen, M the mass of the 
sample, m is the mass of the mobile crew. 


Phase angle, after a rapid increasing, looks having a constant 
behavior during the test. Almost for the whole test this 
values is about 47 degrees. (See Fig. 8). 
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Figure 7 phase angle evolution during a fatigue test (160pe) 


a) Two Point bending test Procedure 

Same loading conditions were considered for 2 Point bending 
test in Nantes. Fatigue tests were undertaken in controlled 
displacement mode at 15 Hz and 20 °C (en- 12697-24-2007, 
EN-12697-26-2007)[16]. 

Strain levels between 120-190 pe were chosen for the 2PB 
tests. Typical strain wave form applied to the trapezoidal 
specimen and typical stress wave form recorded are in Fig 6a 
and Fig 6b. 



(a) (b) 

Figure 5 (a) Strain waveform (b) Stress waveform 


Fig. 9 shows the decreasing of the dissipated energy for 
controlled displacement mode undertaken with the 2PB; 
instead Figure 10 shows the evolution of the hysteresis loop 
(dissipated energy) during a fatigue test. 
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Figure 8 Dissipated Energy versus Number of cycle during a 
fatigue test (170pe) fitted with power law 
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Figure 9 Evolution of hysteresis loop during a fatigue test 



III. Conclusion 

Fatigue tests are undertaken on the same material, at the same 
loading conditions, same frequency (15 Hz) and same 
temperature (20°C); however differences are observed in the 
determination of stiffness, phase angle, number of cycles of 
failure and dissipated energy. 

Several reasons may explain those behaviours : 

- in 2 PB and shear stresses exist in the specimen; 

- different boundary conditions, fixed trapezoidal specimen 

that is glued at the top and at the bottom between two 
plates in 2PB ; 

- presence of shear forces in the 2PB bending test; 

- other phenomena such as heating can influence the result 

differently in the two kind of tests 

The most important point is to understand which testing 
machine simulated better the real fatigue phenomenon; which 
machine gives us the right response of the material, those are 
not easy questions, thus further analysis within a wider 
spectrum of temperatures and frequencies to them is needed. 
Furthermore, dissipated energy ( depending on strain stress 
and phase angle) is a good parameter to describe fatigue 
behavior of bituminous materials and the key point to better 
understand the mechanical response of the same material for 
different kinds of the test. 

1 ) the complex modulus test 


a) introduction 

The test module allows the complex to model empirically the 
viscoelastic behavior of the asphalt-coated under cyclic 
loading. 

The specimens are subjected to repeated stresses under 
zero-centered sinusoidal load. The sample complex modulus 
is measured at different frequencies and at different 
temperatures. These two parameters are set for each 
individual measure. Frequencies and temperatures usually 
tested depend on the test apparatus and capabilities of the 
materials tested; they vary between -20 ° C and 40 ° C and that 
of the frequency from 1 to 40 Hz. 

According to the NF-EN- 12697-26 standard, the 
deformations applied to asphalt specimen must be kept below 
50, 10 6 m / m to prevent fatigue damage to the mix specimen 
tested. Note that by increasing the number of load cycles (up 
to several thousands) on the same complex modulus test 


(homogeneous or inhomogeneous) it leads to another type of 
test: fatigue test. 

b ) Experim ental study 

E L is a complex number composed of a real part E x and an 
imaginary part E 2 . The E L can also be divided into two 
parameters in its vector form; its amplitude ( E L ) and 
argument (c|>). 



Figure 10: Representation of the complex modulus 


The representations of the most used results of the complex 
modulus are those in the complex plane E 2 = ftE^i also 
called representation Cole and Cole (Fig. 12), In the space of 
Black S = fi[|E"D called Black curve (Fig. 13) and 
representation of the main curve. E = f(a T . fj (Fig.l4)ata 
given reference temperature. 



Figure 11: Representation in the Cole and Cole plan 


B. 1.4.1 Measuring criteria 

For asphalt mixes, horizontal bottom layers’ distortion is 
written by a law of the type: 

^ = ^ 6 (Af" L '(20) 

with — V), of 5 th order 



Dynamic Modulus |£*( (MPa) 

Figure 12: Black representation in the plan 
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1 ) representations of complex modulus results 

a) Analysis of results in the complex plane and in the 
space Black 

In Fig. 13, which represents the variation of the viscous 
modulus as a function of the elastic modulus, it appears that 
the experimental points are placed very substantially on the 
same global curve, independent of temperature and 
frequency. On this curve, only one dependent on these factors 
the curvilinear abscissa of the experimental point. The 
important consequence demonstrated by Huet[ 12] of this 
representation is that the frequency and temperature play 
interchangeable roles with respect to complex modulus. 

This curve tends toward the origin for increasing temperatures 
and decreasing frequencies, and towards an asymptotic point 
E s located on the real axis when the frequency approaches 
infinity and when the temperature decreases. Materials 
behave increasingly as pure elastic materials when the 
frequency becomes very great or when the temperature 
becomes very low. The shape of the curve in Fig. 13 is roughly 
that of a circular arc, passing by E-. (near 0) and E^and 
having a significant asymmetry. Originally, the angle that the 
experimental curve makes with E x axis is much greater than 
what it makes with this axis at the other end. The fact that the 
curve seems to pass through the origin indicates a liquid 
behavior. This behavior has been observed on creep tests. 
This curve can be used to calibrate the rheological models. 
The representation in the complex plane, if it is accurate, 
however, has a drawback. Indeed, if coordinates are 
arithmetic, accuracy becomes lower for low values of the 
module. To improve calibration at low modulus values, we 
use most often the representation in the Black space. In this 
new representation, the experimental points fall around the 
same overall curve. The extrapolation of the curve to E-, is 
much easier than the complex plane. This curve is used to 
display the low modulus area and also to reduce of the phase 
angle for the high temperatures. The argument of the complex 
modulus also known as phase angle believes under the same 
conditions, from 0 to a finite value below tt/ 2. At constant 
temperature, the reverse occurs when the frequency increases. 
In this case the material is considered as Theologically simple. 


h ) Master curve Analysis 

When analyzing the isothermal curves shown in Fig. 14, we 
see that the same value of E L can be obtained for different 
pairs (frequency, temperature). This property is called 
time-temperature equivalence property [17]. It is translated 
by writing E “(te, Tl in the form E"Ctdf(T']j . When the 
principle of time-temperature equivalence has been observed 
for a given material, it is possible to construct a single curve 
(logEMcgfD , where F is the frequency of an arbitrarily 
chosen reference temperature T R . 

This curve is obtained by translation parallel to the frequency 
axis of each isotherm compared with the isotherm 
corresponding to the reference temperature till the 
superposition of the same order points. 

Thus the obtained curve is called the master curve at 
reference temperature. It allows obtaining values of modulus 
at frequencies inaccessible by experimentation. It is generally 


considered a reference temperature 15°C for coated-materials 
for the construction of this master curve. 

100020 


ioooo 

Cl 

51000 

u 

100 


10 

0.0CO1 0.01 1 100 IOOOO 1000000 

Equivalent Frequency (Hz) 

Figure 13: Design principle of a master curve at 15°C 



c) Determination of the translation factor 

The equivalence relation (frequency-temperature) was 
observed experimentally long before a theoretical 
justification had been given for this phenomenon, in fact, It 
was found experimentally that the effects of temperature and 
frequency are inverse to one another and it is possible to 
obtain in most cases the same results values by choosing 
properly the values of frequency and temperature. 

This occurs on the experimental results obtained at different 
temperatures between the glass transition point and the 
rubbery transition point in that it is possible to bring these 
curves to each other by means of a translation. It is possible to 
empirically determine the values of these translations as a 
function of temperature and thereby know the equivalence 
relation or the translation factor relative to a given reference 
temperature a T . The translation factor makes it possible to 
move from a test at a temperature and frequency f L to a test 
conducted at temperature T 0 and a frequency f L =£ f and equal 
to . Various forms for the expression of have been 
proposed. The most used at present is the relationship of 
William-Landel-F erry (WLF) [17]. 

This relationship has been proposed from the experimental 
results obtained on different polymers, before a theoretical 
justification based on molecular chains theory is attributed by 
Rouse[18]. 

The translation coefficient of isotherm T compared with the 
isotherm of chosen reference is log(a T 3 such that 
E"(w,T 3 = E"(ay. a-pTftj with i Tr = 1. 

Several formulas are available to determine log a T . The most 
used are: 

• the Arrhenius equation: 

= ?£-£) (2D 

AH is the apparent activation energy depending on the 
material; R is the constant of perfect gas; T and are the 
temperatures expressed in Kelvin, 

• the relationship of William-Landel-F erry (WLF): 

from the analysis of the rheological results on 
different materials [19] they proposed a law of 

translation factors variation log(a T ) between the 
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experimental isotherms T and a temperature of any 


reference 


log(a T } = 


-CaCT-Tj.) 

Ci+CT-Tjf) 


( 22 ) 


and C z parameters are constants depending on the material 
studied and have the following properties: 

Cl- Cj = c; . C 2 ; (T raf - Cj > = CT„ f - q > (23) 

T and TR are the temperatures expressed in Kelvin; and C z 
are obtained at T r€f while and C z are obtained at T r ' € f[19]. 
An example of changing the translation coefficient obtained 
from the WLF relationship depending on the temperature is 
shown in Fig. 15 



Figure 14: translation coefficient as a function of temperature 
T rs f = 1 5 : C with C x = 32 and C z = 216 of the WLF law 
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• The translation factor can also be described by a 
second order polynomial 

l dq (*St 1 = aT 2 + b'T + c (24) 
a, b, c are regression coefficients 
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